Efficient Anchoring of Erianthus arundinaceus Chromatin Introgressed into Sugarcane by Specific Molecular Markers

Erianthus arundinaceus is a valuable gene reservoir for sugarcane improvement. However, insufficient molecular markers for high-accuracy identification and tracking of the introgression status of E. arundinaceus chromatin impede sugarcane breeding. Fortunately, suppression subtractive hybridization (SSH) technology provides an excellent opportunity for the development of high-throughput E. arundinaceus-specific molecular markers at a reasonable cost. In this study, we constructed a SSH library of E. arundinaceus. In total, 288 clones of E. arundinaceus-specific repetitive sequences were screened out and their distribution patterns on chromosomes were characterized by fluorescence in situ hybridization (FISH). A subtelomeric repetitive sequence Ea086 and a diffusive repetitive sequence Ea009, plus 45S rDNA-bearing E. arundinaceus chromosome repetitive sequence EaITS were developed as E. arundinaceus-specific molecular markers, namely, Ea086-128, Ea009-257, and EaITS-278, covering all the E. arundinaceus chromosomes for high-accuracy identification of putative progeny. Both Ea086-128 and Ea009-257 were successfully applied to identify the authenticity of F1, BC1, BC2, BC3, and BC4 progeny between sugarcane and E. arundinaceus. In addition, EaITS-278 was a 45S rDNA-bearing E. arundinaceus chromosome-specific molecular marker for rapid tracking of the inherited status of this chromosome in a sugarcane background. Three BC3 progeny had apparently lost the 45S rDNA-bearing E. arundinaceus chromosome. We reported herein a highly effective and reliable SSH-based technology for discovery of high-throughput E. arundinaceus-specific sequences bearing high potential as molecular markers. Given its reliability and savings in time and efforts, the method is also suitable for development of species-specific molecular markers for other important wild relatives to accelerate introgression of wild relatives into sugarcane.

Development of SLAF-Sequence and Multiplex SNaPshot Panels for Population Genetic Diversity Analysis and Construction of DNA Fingerprints for Sugarcane.

A genetic diversity analysis and identification of plant germplasms and varieties are important and necessary for plant breeding. Deoxyribonucleotide (DNA) fingerprints based on genomic molecular markers play an important role in accurate germplasm identification. In this study, Specific-Locus Amplified Fragment Sequencing (SLAF-seq) was conducted for a sugarcane population with 103 cultivated and wild accessions. In total, 105,325 genomic single nucleotide polymorphisms (SNPs) were called successfully to analyze population components and genetic diversity. The genetic diversity of the population was complex and clustered into two major subpopulations. A principal component analysis (PCA) showed that these accessions could not be completely classified based on geographical origin. After filtration, screening, and comparison, 192 uniformly-distributed SNP loci were selected for the 32 chromosomes of sugarcane. An SNP complex genotyping detection system was established using the SNaPshot typing method and used for the precise genotyping and identification of 180 sugarcane germplasm samples. According to the stability and polymorphism of the SNPs, 32 high-quality SNP markers were obtained and successfully used to construct the first SNP fingerprinting and quick response codes (QR codes) for sugarcane. The results provide new insights for genotyping, classifying, and identifying germplasm and resources for sugarcane breeding.

Gene Co-Expression Analysis Reveals Transcriptome Divergence between Wild and Cultivated Sugarcane under Drought Stress.

Drought is the main abiotic stress that constrains sugarcane growth and production. To understand the molecular mechanisms that govern drought stress, we performed a comprehensive comparative analysis of physiological changes and transcriptome dynamics related to drought stress of highly drought-resistant (ROC22, cultivated genotype) and weakly drought-resistant (Badila, wild genotype) sugarcane, in a time-course experiment (0 h, 4 h, 8 h, 16 h and 32 h). Physiological examination reviewed that ROC22, which shows superior drought tolerance relative to Badila, has high performance photosynthesis and better anti-oxidation defenses under drought conditions. The time series dataset enabled the identification of important hubs and connections of gene expression networks. We identified 36,956 differentially expressed genes (DEGs) in response to drought stress. Of these, 15,871 DEGs were shared by the two genotypes, and 16,662 and 4423 DEGs were unique to ROC22 and Badila, respectively. Abscisic acid (ABA)-activated signaling pathway, response to water deprivation, response to salt stress and photosynthesis-related processes showed significant enrichment in the two genotypes under drought stress. At 4 h of drought stress, ROC22 had earlier stress signal transduction and specific up-regulation of the processes response to ABA, L-proline biosynthesis and MAPK signaling pathway-plant than Badila. WGCNA analysis used to compile a gene regulatory network for ROC22 and Badila leaves exposed to drought stress revealed important candidate genes, including several classical transcription factors: NAC87, JAMYB, bHLH84, NAC21/22, HOX24 and MYB102, which are related to some antioxidants and trehalose, and other genes. These results provide new insights and resources for future research and cultivation of drought-tolerant sugarcane varieties.

Chromosome-specific painting unveils chromosomal fusions and distinct allopolyploid species in the Saccharum complex.

Karyotypes provide key cytogenetic information on the phylogenetic relationships and evolutionary origins in related eukaryotic species. Despite our knowledge of the chromosome numbers of sugarcane and its wild relatives, the chromosome composition and evolution among the species in the Saccharum complex have been elusive owing to the complex polyploidy and the large numbers of chromosomes of these species. Oligonucleotide-based chromosome painting has become a powerful tool of cytogenetic studies especially for plant species with large numbers of chromosomes. We developed oligo-based chromosome painting probes for all 10 chromosomes in Saccharum officinarum (2n = 8x = 80). The 10 painting probes generated robust fluorescence in situ hybridization signals in all plant species within the Saccharum complex, including species in the genera Saccharum, Miscanthus, Narenga and Erianthus. We conducted comparative chromosome analysis using the same set of probes among species from four different genera within the Saccharum complex. Excitingly, we discovered several novel cytotypes and chromosome rearrangements in these species. We discovered that fusion from two different chromosomes is a common type of chromosome rearrangement associated with the species in the Saccharum complex. Such fusion events changed the basic chromosome number and resulted in distinct allopolyploids in the Saccharum complex.

Genome-wide identification and expression analysis of AP2/ERF transcription factors in sugarcane (Saccharum spontaneum L.).

BACKGROUND: APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) transcription factors play essential roles in plant growth, development, metabolism, and responses to biotic and abiotic stresses. However, few studies concerning AP2/ERF genes in sugarcane which are the most critical sugar and energy crops worldwide. RESULTS: A total of 218 AP2/ERF genes were identified in the Saccharum spontaneum genome. Phylogenetic analysis showed that these genes could be divided into four groups, including 43 AP2s, 160 ERFs and Dehydration-responsive element-binding (DREB) factors, 11 ABI3/VPs (RAV), and four Soloist genes. These genes were unevenly distributed on 32 chromosomes. The structural analysis of SsAP2/ERF genes showed that 91 SsAP2/ERFs lacked introns. Sugarcane and sorghum had a collinear relationship between 168 SsAP2/ERF genes and sorghum AP2/ERF genes that reflected their similarity. Multiple cis-regulatory elements (CREs) present in the SsAP2/ERF promoter were related to abiotic stresses, suggesting that SsAP2/ERF activity could contribute to sugarcane adaptation to environmental changes. The tissue-specific analysis showed spatiotemporal expression of SsAP2/ERF in the stems and leaves of sugarcane at different development stages. In ten sugarcane samples, 39 SsAP2/ERFs were not expressed, whereas 58 SsAP2/ERFs were expressed in all samples. Quantitative PCR experiments showed that SsERF52 expression was up-regulated under salt stress, but suppressed under dehydration stress. SsSoloist4 had the most considerable upregulation in response to treatment with the exogenous hormones ABA and GA. Within 3 h of ABA or PEG6000 treatment, SsSoloist4 expression was up-regulated, indicating that this gene could play a role in the responses to ABA and GA-associated dehydration stress. Analysis of AP2/ERF gene expression patterns under different treatments indicated that SsAP2/ERF genes played an essential role in dehydration and salt stress responses of S. spontaneum. CONCLUSIONS: In this study, a total of 218 members of the AP2 / ERF superfamily were identified in sugarcane, and their genetic structure, evolution characteristics, and expression patterns were studied and analyzed. The results of this study provide a foundation for future analyses to elucidate the importance of AP2/ERF transcription factors in the function and molecular breeding of sugarcane.

Sequence Evolution, Abundance, and Chromosomal Distribution of Ty1-copia Retrotransposons in the Saccharum spontaneum Genome.

Saccharum spontaneum is a wild germplasm resource of the genus Saccharum that has many valuable traits. Ty1-copia retrotransposons constitute a large proportion of plant genomes and affect genome sequence organization and evolution. This study aims to analyze the sequence heterogeneity, phylogenetic diversity, copy number, and chromosomal dispersion patterns of Ty1-copia retrotransposons in S. spontaneum. A total of 44 Ty1-copia reverse transcriptase subclones isolated from S. spontaneum showed a range of heterogeneity, and all sequences were A-T rich, averaging approximately 54.59%. Phylogenetic analysis divided the 44 reverse transcriptase sequences into 5 distinct lineages (Retrofit/Ale, Sire/Maximus, Bianca, Tork/TAR, and Ty1-copia like). Dot-blot hybridization revealed that Ty1-copia retrotransposons consisted of a significant component of approximately 38,900 copies and 16,300 copies per genome in the accessions YN82-114 (2n = 10x = 80) and AP85-441 (2n = 4x = 32), respectively. The results of a local blast analysis showed that there are 15,069 Ty1-copia retrotransposon copies in the genome of AP85-441, of which the Retrofit/Ale lineage had the highest copy number, followed by the Tork/TAR, Sire/Maximus, and Bianca lineages. Furthermore, both FISH and the local blast analysis with AP85-441 genomic data demonstrated that the Ty1-copia retrotransposons were unevenly distributed throughout the chromosomes. Taken together, this study provides insights into the role of Ty1-copia retrotransposons in the evolution and organization of the S. spontaneum genome.

Breeding signature of combining ability improvement revealed by a genomic variation map from recurrent selection population in Brassica napus.

Combining ability is crucial for parent selection in crop hybrid breeding. The present investigation and results had revealed the underlying genetic factors which might contribute in adequate combining ability, further assisting in enhancing heterosis and stability. Here, we conducted a large-scale analysis of genomic variation in order to define genomic regions affecting the combining ability in recurrent selection population of rapeseed. A population of 175 individuals was genotyped with the Brassica60K SNP chip. 525 hybrids were assembled with three different testers and used to evaluate the general combining ability (GCA) in three environments. By detecting the changes of the genomic variation, we identified 376 potential genome regions, spanning 3.03% of rapeseed genome which provided QTL-level resolution on potentially selected variants. More than 96% of these regions were located in the C subgenome, indicating that C subgenome had sustained stronger selection pressure in the breeding program than the A subgenome. In addition, a high level of linkage disequilibrium in rapeseed genome was detected, suggesting that marker-assisted selection for the population improvement might be easily implemented. This study outlines the evidence for high GCA on a genomic level and provided underlying molecular mechanism for recurrent selection improvement in B. napus.

Genetic dissection of heterosis using epistatic association mapping in a partial NCII mating design.

Heterosis refers to the phenomenon in which an F1 hybrid exhibits enhanced growth or agronomic performance. However, previous theoretical studies on heterosis have been based on bi-parental segregating populations instead of F1 hybrids. To understand the genetic basis of heterosis, here we used a subset of F1 hybrids, named a partial North Carolina II design, to perform association mapping for dependent variables: original trait value, general combining ability (GCA), specific combining ability (SCA) and mid-parental heterosis (MPH). Our models jointly fitted all the additive, dominance and epistatic effects. The analyses resulted in several important findings: 1) Main components are additive and additive-by-additive effects for GCA and dominance-related effects for SCA and MPH, and additive-by-dominant effect for MPH was partly identified as additive effect; 2) the ranking of factors affecting heterosis was dominance > dominance-by-dominance > over-dominance > complete dominance; and 3) increasing the proportion of F1 hybrids in the population could significantly increase the power to detect dominance-related effects, and slightly reduce the power to detect additive and additive-by-additive effects. Analyses of cotton and rapeseed datasets showed that more additive-by-additive QTL were detected from GCA than from trait phenotype, and fewer QTL were from MPH than from other dependent variables.

Interacted QTL mapping in partial NCII design provides evidences for breeding by design.

The utilization of heterosis in rice, maize and rapeseed has revolutionized crop production. Although elite hybrid cultivars are mainly derived from the F1 crosses between two groups of parents, named NCII mating design, little has been known about the methodology of how interacted effects influence quantitative trait performance in the population. To bridge genetic analysis with hybrid breeding, here we integrated an interacted QTL mapping approach with breeding by design in partial NCII mating design. All the potential main and interacted effects were included in one full model. If the number of the effects is huge, bulked segregant analysis were used to test which effects were associated with the trait. All the selected effects were further shrunk by empirical Bayesian, so significant effects could be identified. A series of Monte Carlo simulations was performed to validate the new method. Furthermore, all the significant effects were used to calculate genotypic values of all the missing F1 hybrids, and all these F1 phenotypic or genotypic values were used to predict elite parents and parental combinations. Finally, the new method was adopted to dissect the genetic foundation of oil content in 441 rapeseed parents and 284 F1 hybrids. As a result, 8 main-effect QTL and 37 interacted QTL were found and used to predict 10 elite restorer lines, 10 elite sterile lines and 10 elite parental crosses. Similar results across various methods and in previous studies and a high correlation coefficient (0.76) between the predicted and observed phenotypes validated the proposed method in this study.